import math as m
import turtle as t


#计算太阳赤纬



t.penup()



h2=1
a3=2
t.goto(a3,h2)
t.pendown()
t.done()
mm=input("1请输入天数")
n=int(mm)  #把输入的字符转化为数据

for i in range(1,13,1):
    alpha=0.006918-0.399912*m.cos((2*m.pi*(n-1))/365)+0.070257*m.sin((2*m.pi*(n-1))/365)-0.006758*m.cos(2*((2*m.pi*(n-1))/365))+0.000907*m.sin(2*((2*m.pi*(n-1))/365))-0.002697*m.cos(3*((2*m.pi*(n-1))/365))+0.00148*m.sin(3*((2*m.pi*(n-1))/365))  #用公式计算
    beta=alpha/(0.0174533)  #把计算结果从弧度转化为角度
    print("当日的太阳赤纬值为：")
    print(beta)

 
    #录入北京的纬度信息并转化为角度
    y=((m.pi)/180)*(39+(54/60)+(20/3600)) 


    #计算太阳时角
    tt=input("请输入时间")
    t1=float(tt)  #把输入的字符转化为数据
    print("请查时差表")
    ttt=input("时差（分钟）")
    tttt=input("时差（秒钟）")
    mmm=int(ttt)  #把输入的字符转化为数据
    s=int(tttt)  #把输入的字符转化为数据
    x=mmm/60+s/3600   #把角度转化为100进制
    t2=t1-(14+(44/60))/60  #把角度转化为100进制
    t3=t2+x
    t=(((m.pi)/180)*15)*(t3-12)   #代入公式


    #求太阳的高度角
    gamma=m.sin(y)*m.sin(alpha)+m.cos(y)*m.cos(alpha)*m.cos(t)  #代公式求值
    hhh=m.asin(gamma)  #同上
    h=hhh/(0.0174533)  #把计算结果从弧度转化为角度
    print("当时该时太阳的高度角为：")
    print(h)


    #求太阳的方位角
    a1=(m.sin(alpha)-m.sin(hhh)*m.sin(y))/(m.cos(hhh)*m.cos(y))
    a2=m.acos(a1)
    a=a2/(0.0174533)
    print("当时太阳的方位角为：")
    print(a)


